Top anti-top pairs at the LHC heavy ion collision: a new interesting probe of quark 

gluon plasma 
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We investigate top anti-top quark pair production in lead-lead collisions at the Large Hadron 
Collider with nucleon-nucleon center of mass energy of 5.5 TeV. Due to the very high temperature 
and energy density created in heavy ion collision, a new state of QCD matter known as Quark- 
Gluon Plasma (QGP) is expected to be produced. Top decay products loose energy inside the QGP 
medium. Therefore, we also study the medium modifications of different kinematic distributions. 
We observe significant modification in the dijets and trijets invariant mass distributions. We also 
found that the peak position and shape of the distributions could be used to characterize the nature 
of jet energy loss in the QGP. 

PACS numbers: 



The primary goal of heavy ion collisions (HIC) at Rel- 
ativistic Heavy Ion Collider (RHIC) at BNL, and the 
Large Hadron Collider (LHC) at CERN is to produce 
and study the properties of a hot/dense state of QCD 
matter known as Quark-Gluon Plasma (QGP) QGP 
is a deconfined state of matter where quarks and gluons 
are the effective degrees of freedom rather than nucleons 
or hadrons Fast partons propagating in a hot/dense 
nuclear medium are expected to loose a large fraction 
of their energy Q. The observation of the suppression 
of energetic partons in the QGP, that is jet quenching 

, and centrality-dependent dijet asymmetry fsl are the 
most important results from the HIC at RHIC and LHC 
experiment, respectively. 

M^/Z-bosons are massive weakly interacting Standard 
Model (SM) particles. The larger LHC HIC energies 
open the possibility to probe the nucleus-nucleus colli- 
sions via the PF/Z-bosons. The vector bosons are pro- 
duced early (l/Mw^z) ^ 10~^ fm/c) and their decay 
time is small (tz ^ 0.08 fm/c and tw ^ 0.09 fm/c 0). 
Whereas, in the most accepted picture of the QGP for- 
mation and evolution, at the LHC, QGP is expected to 
form after 1/Aqcd ~ 1 fm/c of the initial hard scat- 
tering, thermalize quickly and it might last ~ 10 fm/c. 
Therefore, weak bosons are produced and decay before 
the formation of QGP and the decay products of pass 
through the QGP. W/Z-hoaon dominantly decays to a 
pair of quarks. Quarks loose energy in the QGP and thus 
the hadronic decays of W/Z-hoson could be an interest- 
ing probe to characterize QGP. However, in presence of 
huge QCD dijet background, it is extremely challenging 
to study the hadronic decays of W^/Z-bosons. The lep- 
tonic BF of W/Z-hoson is small however, due to small 
background, the signature could be easily detected at the 
LHC HIC. Unlike jets, leptons interact electromagneti- 
cally with the QGP and loose experimentally insignifi- 



cant amount of energy within the QGP (TJ. Therefore, 
the leptonic decays of weak bosons behave as a medium 
blind reference. ATLAS collaboration has already mea- 
sured the W/Z-hoson yield in the leptonic decay channels 
for ^/IWn =2.76 TeV @. 

A new regime of heavy ion physics will be reached at 
the LHC with ^/snn = 5.5 TeV where hard and semi- 
hard particle production can dominate over the under- 
lying soft events. The higher LHC energies open the 
possibility to study top quarks for the first time at the 
HIC. The top quark was discovered at the Tevatron ex- 
periment . After the discovery, properties of top quark 
have been extensively studied at the Tevatron and LHC. 
Presently, different properties of top quark is known with 
good precision. As an example, top quark mass and full 
decay width are determined to be mt = 173.5 ± 0.6 GeV 
and Ft = 2.0lg'g GeV Q. Different decay channels and 
branching ratios of top quark have also been observed 
and measured. As a result, top quark can now be con- 
sidered as standard benchmark for other experimental 
observations. In view of this fact, it is important to 
investigate top quarks at the LHC HIC and study the 
medium (QGP) influence on different kinematic distri- 
butions which are precisely known from the previous pp 
and pp collider experiments. In this letter, we have for 
the first time studied ti production at the LHC HIC and 
proposed few kinematic distributions for the study of the 
QGP created in the LHC HIC. 

The main source of top quarks at the LHC HIC is the 
top anti-top (tt) pair production. At leading order in per- 
turbation theory there are two processes that contribute 
to tt production: quark- antiquark annihilation, qq tt 
and gluon-gluon fusion, gg — ^ ti. With 5.5 TeV center-of- 
mass energy per nucleon, the NLO-I-NNLL ti production 
cross-section per nucleon-nucleon collision in Pb-Pb re- 
action is given hy aNN (ti) = 80.6 pb 0. Therefore, the 
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total tt production cross-section in minimum bias Pb-Pb 
scattering is estimated to be 3.5 /x6 in the frame work of 



the Glauber model The instantaneous luminosity 



of the LHC Pb-Pb collision at -v/s^vAf = 5.5 TeV is ex- 
pected to be lO^^cm^s"^. Therefore, 1 nb""'^ integrated 
luminosity data will be accumulated within one month 
(10^ second) of Pb-Pb colhsion at ^/sJ^J^ = 5.5 TeV. 
These estimations indicate that about 3500 tt events are 
expected to be produced in one month of LHC HIC run- 
ning at v^^jv = 5.5 TeV. 

Due to the large top mass, tt pairs are produced early 
in the HIC. The decay width of top quark is large and 
thus the decay time is small. As a result, top quarks are 
produced and decay before formation of QGP. Top quark 
decays to a bottom quark (&) and VF-boson with almost 
100% branching fraction: t — bW^ . VF-boson subse- 
quently decays hadronically {W^ qq') with 67.7% BE 
or leptonically {W^ In) with 32.3% BF &]. Therefore, 
for tt production, there are only three possible final state 
topologies: (i) Hadronically decaying tt pairs: Both the 
top quarks decay hadronically (t — bqq') and give rise to 
2-6 jets and four light quarks jets in the final state with 
46% effective branching ratio, (ii) Semi-leptonically de- 
caying tt pairs: One top quark decays hadronically and 
the other decays leptonically (i — >■ bliyi). The final state is 
characterized by 2-6 jets-|-2-light quark jets-|-one charged 
lepton -I- one neutrino. The effective branching fraction 
of semi-leptonic decay mode of tt pairs is 43.7%. (Hi) 
Leptonically decaying tt pairs: In this case, both the top 
quarks decay leptonically giving rise to 2-6 jets -I- 2-lepton 
-I- 2-neutrino final state with 10.4% effective branching 
fraction. 

The jets (6-jets as well as light jets) and leptons can be 
observed at the LHC. The neutrinos remain invisible at 
the detectors and give rise to a imbalance in the trans- 
verse momentum known as missing transverse momen- 
tum. However, in the HIC environment, faithful mea- 
surement of missing transverse momentum will be chal- 
lenging in the presence of a continuum energy deposit 
from the produced QGP. 

Hadronically decaying ti events suffer from huge QCD 
background. Moreover, due to a large combinatorial 
background it is difficult to reconstruct top quark and 
ly-boson from hadronically decaying tt pairs. Lepton- 
ically decaying tt channel is a clean signal channel due 
to less background. However, the rate of this channel 
is suppressed by the top quark leptonic branching frac- 
tion. Moreover, in the absence of proper knowledge about 
missing transverse momentum, the reconstruction of top 
quark and M^-boson mass will be difficult from tt lep- 
tonic decay channel. As a result, in this letter, we have 
investigated the semi-leptonic decay channel of tt pairs 
due to the following advantages: (i) The rate of semi- 
leptonic tt final state is relatively large. If we consider 
only electron and muon decay modes, Pb-Pb collision 
with -y/s^Tv = 5.5 TeV gives rise to about 1000 semi- 



leptonically decaying tt events for 1 nb""'^ integrated lu- 
minosity, (ii) Due to the presence of a lepton, the semi- 
leptonic final state of tt pairs suffers less from the QCD 
background, (iii) At the parton level, the semi-leptonic 
tt final state contains only 2-light quark jets arising from 
the VF-decay. Therefore, W invariant mass can be con- 
structed with out any ambiguity. However, there is a two 
fold ambiguity in the reconstruction of top quark mass. 

After the decay of tt pairs, the decay products pass 
through the hot and dense QGP medium and thus, loose 
energy. The energy loss of energetic partons, so called 
"jet quenching" , leads to a number of phenomena which 
are already seen at RHIC and LHC. In this work, we 
have investigated the quenching of ti decay products and 
proposed some new phenomena which could be seen at 
the LHC HIC. We have used PYTHIA [H to simu- 
late the production and decay of tt pairs. Subsequently, 
the PYTHIA generated tt events are passed in to a fast 
Monte-Carlo simulator PYQUEN [l3| for simulating the 
energy loss (quenching) of tt decay products. Finally, 
quarks and gluons are hadronized according to the Lund 
string mode PYQUEN simulates the radiative and coUi- 
sional (Coll.) energy loss [1^ of hard partons in longitu- 
dinally expanding QGP taking into account the realistic 
nuclear geometry. The radiative energy loss is calculated 
in the framework of BDMS model [l5| with the simple 
generalization to a massive quark case using the dead- 
cone approximation TGj. Measuring jet energy as a sum 
of the energies of final hadrons moving inside an angu- 
lar cone with a given finite size allow some of the radi- 
ated gluons to belong to the jet and thus some part of 
the radiated energy to be reconstructed. Therefore, the 
knowledge of angular structure of medium-induced radi- 
ation is very important for any phenomenological predic- 
tion using jets. In this analysis, we have used the simple 
parametrizations of the gluon distribution over the emis- 
sion angle 6 available in PYQUEN: (i) Small-angular ra- 
diation (SAR): dN9/de cx sin6lexp[-(6'-6'o)V26'o]> where 
9o ^ 5° is the typical angle of the coherent gluon radia- 
tion as estimated in Ref. [17| ; (ii) Wide-angular radiation 
(WAR): dNS/de cx 1/6. The strength of the energy loss 
in PYQUEN is determined mainly by the initial maxi- 
mal temperature T™°^ of hot matter in Pb-Pb collisions. 
The energy loss also depends on the number N f oi 'Ac- 
tive fiavors in the medium. In our analysis, we have used 

j^vaax ^ ^ Nf =2. 

In this letter, we have investigated 4-jets out of which 
2-jets are 6-tagged plus one lepton signature as signal of 
tt production in HIC. Therefore, before going into the 
details of our analysis, it is important to discuss the sta- 
tus of jet reconstruction and 6-tagging in the context of 
HIC. The main obstacle to studying jets in HIC is the 
presence of the huge background given by the underly- 
ing event (UE) This UE needs to be properly subtracted 
from the momentum of a given jet in order to recon- 
struct its true momentum. It was shown in Ref. [l8| that 



3 



> 

a) 
P 




"T 1 1 1 1 — I — r 

Coll.+WAR 
C0II.+ SAR 



Background x 5'^ 



85 90 50 55 60 65 70 75 80 85 90 



mjj [GeV] 



FIG. 1: Invariant mass distribution of hardest and second 
hardest non &-jets for different energy loss scenarios. To de- 
termine the peak position of the distribution, we have fitted 
these distributions with asymmetric Gaussian functions. The 
parameters of fitting are presented in Table. |I] 
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FIG. 2: Same as Fig. [T] but for nibjj distribution. 



in presence of QGP, faithful reconstruction of jets are 
possible using anti-A:j- algorithm. Jets have already been 
successfully used as an observable in HIC by the AT- 
LAS collaboration Moreover, in Ref. [19|, different 
jet shape variables have been studied for the jets passing 
through QGP. In Ref ^2^, 6-tagging in the environment 
of HIC has been studied by examining the reconstruction 
efficiency and rejection power (against light quark jets) 
using secondary vertex finding. Their study suggests that 
Cb — 50% 6-tagging efficiency can be achieved at the HIC 
for a rejection power of 50. 

In our analysis, we have introduced a set of basic 
selection criteria to identify electrons, muons, jets etc. 
The object selection is described in brief in the follow- 
ing: (i) Jets are constructed using anti-fcy algorithm with 
R = 0.4 and only jets with pr > 20 GeV and I77I < 2.5 
are considered for further analysis. To take into account 
the effects of finite detector resolution, we have smeared 
the jet energies with Gaussian functions in Ref. 2l|. (ii) 



TABLE I: Parameters used for the Gaussian fitting of 
Fig [T] and [1 
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We demand that Icpton candidates (both electron and 
muon) have pr > 20 GeV and are separated from jets 
by at least AR — 0.5. After reconstructing different 
objects, we consider events with one lepton (electron or 
muon) and > 4-jets for further analysis. We also demand 
that out of the 4-jets, two jets are ^-tagged. The dom- 
inant background for semileptonic tt signal arises from 
W/Z+jets production followed by the leptonic decay of 
W/Z-hoson. Here, &-jets results from the mistagging of 
light jets. Since the mistagging efficiency of light quark 
jets to be tagged as 6-jets is small, 5-tagging significantly 
reduce this background. Production of W/Z 55-f jets also 
contributes to the background. However, we have esti- 
mated that these cross-sections are very small compared 
to the tt cross-section. Therefore, semileptonic decay 
products of tt could be easily detected over the back- 
ground. 

To study the infiuence of the hot/dense QGP, we have 
constructed the following kinematic distributions. 
Di-jets invariant mass distribution: We have ordered the 
non 6 jets according to their pT hardness (p^^ > Pt > ■•■) 
and constructed the invariant mass of the hardest (ji) 
and the second hardest (J2) jet, mjj. In the semileptonic 
decay of tt pairs, non 5-jets arise from the decay of one 
W-boson. Therefore, in absence of QGP, mjj distri- 
bution should be peaked at the M^-mass, mw = 80.4 
GeV. However, in presence of QGP, VF-decay product 
suffers energy loss. In Fig. [1] we have presented the 
mjj distributions for different energy loss scenarios. To 
determine the peak position of the distributions, we 
have fitted the distributions with asymmetric Gaussian 
functions: 




if m < Too 
otherwise 



where, Oq, TOq and ai_2 are the parameters of fitting. In 
Fig. [T] (left panel), we have presented mjj distributions 
for coUisional and radiative (SAR as well as WAR) energy 
loss scenario separately. Corresponding fitting parame- 
ters are presented in Table HI Fig. [T] (left panel) and Ta- 
ble |T] show that collisional energy loss significantly (about 
5 GeV) shifts the peak position of mjj distribution from 
FK-mass. Whereas, radiative energy loss (both SAR and 
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WAR) gives rise to small change in the peak position of 
nijj distribution. Due to the high boost of hard partons, 
the radiated gluons shift towards the parent parton direc- 
tion and thus, resulting jets include large part of radiated 
gluon energies. As a result, radiative energy loss of hard 
partons has negligible impact on the peak position of 
nijj distribution. However, Fig. [1] (left panel) shows that 
WAR significantly changes the shape (which could be 
quantified by the fitting parameters oq, cti and (T2 in Ta- 
ble |l| of nijj distribution. In the passage of a fast parton 
through QGP, collisional and radiative energy loss occur 
simultaneously. Therefore, in Fig. [1] (right panel), we 
have presented rrijj distribution in presence of both col- 
lisional and radiative energy loss. Corresponding fitting 
parameters are presented in the last two rows of Table HI 
Table |T] shows that the peak position of rrijj distribution 
is determined by the collisional energy loss. Whereas, the 
shape parameters are governed by the nature of radiative 
loss. In Fig. [1] (right panel), we have also presented the 
background. The background contributions are substan- 
tially small compared to the tt contribution. 
Tri-jets (bjj) invariant mass distribution: Hadronic de- 
cay of one top quark gives rise to one &-jet and two light 
quark jets. Therefore, it is important to study invariant 
mass distribution of b jet and di-jets system arising from 
the hadronic top quark decay. However, it is difficult to 
identify the b jet and di-jets arising from the same top 
decay. There are several algorithm available in the lit- 
erature for the reconstruction of top quark. However, 
most of the algorithms rely on the knowledge of ruhjj 
peak position (top mass: rrit = 173.2 GeV). In presence 
of QGP, energy loss of top decay products shift the peak 
position of mbjj distribution. Therefore, most of the top 
reconstruction algorithm are not applicable for Pb-Pb 
collision in their present form. In our analysis, we have 
used the following simplified algorithm for reconstruct- 
ing mbjj peak position. We first order b jets according 
to their pT hardness (p^ > p^) and constructed two 
invariant masses: mb^jj and rnb^jj- Out of these two in- 
variant masses, we consider the invariant mass which is 
close to the nit = 173.2 GeV for plotting mbjj distribu- 
tion. In Fig. [21 we have presented mbjj distributions. In 
Fig. [2] (left panel), we have shown the effect of collisional 
and radiative energy loss on mbjj distribution separately. 
Whereas, right panel of Fig. [2] shows the resulting mbjj 
distribution if we consider both collisional and radiative 
energy loss simultaneously. We have also fitted these dis- 
tributions with asymmetric Gaussian functions and the 
fitting parameters are presented in Table HI Due to the 
energy loss, the position of mbjj peak shifts about 12 
GeV from mj. The shift could be easily observed at the 
LHC as a signature of QGP. 

To summarize, we have investigated semileptonic tt 
signature at the LHC HIC. Semileptonic tt signature 
could be easily observed at the LHC HIC with -\/s^jv ~ 
5.5 TeV. However, due to the presence of QGP tt decay 



products suffer energy loss and thus shape of different 
kinematic distributions are modified significantly. As for 
example, we have studied dijets and trijets invariant mass 
distributions and predicted significant change and shift 
in the shape and peak position of these distributions, re- 
spectively. 
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